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ABSTRACT 

Context. The mass-loss rate is a key parameter of massive stars. Adequate stellar atmosphere models are required for spectral ana- 
lyses and mass-loss determinations. Present models can only account for the inhomogeneity of stellar winds in the approximation 
of small-scale structures that are optically thin. Compared to previous homogeneous models, this treatment of "microclumping" has 
led to reducing empirical mass-loss rates by factors of two to three. Further reductions are presently discussed in the literature, with 
far-reaching consequences e.g. for stellar evolution and stellar yields. 

Aims. Stellar wind clumps can be optically thick in spectral lines. We investigate how this "macroclumping" influences the radiative 
transfer and the emergent line spectra and discuss its impact on empirical mass-loss rates. 

Methods. The Potsdam Wolf-Rayet (PoWR) model atmosphere code is generalized in the "formal integral" to account for clumps 
that are not necessarily optically thin. The stellar wind is characterized by the filling factor of the dense clumps and by their average 
separation. An effective opacity is obtained by adopting a statistical distribution of clumps and applied in the radiative transfer. 
Results. Optically thick clumps reduce the effective opacity. This has a pronounced effect on the emergent spectrum. Our modeling 
for the O-type supergiant f Puppis reveals that the optically thin Hey line is not affected by wind porosity, but that the P v resonance 
doublet becomes significantly weaker when macroclumping is taken into account. The reported discrepancies between resonance-line 
and recombination-line diagnostics can be resolved entirely with the macroclumping modeling without downward revision of the 
mass-loss rate. In the case of Wolf-Rayet stars, we demonstrate for two representative models that stronger lines are typically reduced 
by a factor of two in intensity, while weak lines remain unchanged by porosity effects. 

Conclusions. Mass-loss rates inferred from optically thin emission, such as the Ho- line in O stars, are not influenced by macroclump- 
ing. The strength of optically thick lines, however, is reduced because of the porosity effects. Therefore, neglecting the porosity in 
stellar wind modeling can lead to underestimating empirical mass-loss rates. 
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1. Introduction 

Mass loss by stellar winds plays a key role in the evolution of 
massive stars. The feedback of chemically enriched material, 
momentum, and mechanical energy is important for stellar clus- 
ters and the formation of stars, so it has consequences for the 
whole cosmic evolution (e.g. Leitherer et al. I1992I Woosley et 
al. [1993b . 

It is very important to know the mass-loss rate of stars in their 
different evolutionary phases, and its dependence on initial mass, 
composition, and possibly further parameters. Decades of efforts 
have been needed to establish reliable mass-loss rates for OB 
and Wolf-Rayet (WR) type stars (see Kudritzki & Puis rZOOOb . 
Nevertheless, there are doubts whether the presently accepted 
mass-loss rates are correct. Some recent papers favor a dramatic 
downwards reduction of O-star mass-loss rates with far-reaching 
consequences. The reason for the uncertainties largely stem from 
the implications of wind inhomogeneity, i.e. "clumping". 

One of the numerous proofs of wind clumping is the line 
profile variability routinely found in O and WR spectra, where 
different types of variability can be distinguished. "Discrete ab- 
sorption components" (DACs) are observed regularly in the UV 
spectra from O stars (e.g. Kaper et al. [ 19961). They display pe- 
riodic patterns on timescales of days and are most likely asso- 
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ciated with stellar rotation, but their slow wavelength drift is 
still enigmatic (e.g. Hamann et al. |2001t >. "Modulations", which 
have a slightly shorter timescale and are not strictly periodic, 
are thought to reflect non-radial stellar pulsations (Eversberg et 
al. 11998b . Only the stochastic variability on short timescales of 
hours is attributed to wind clumping. Such stochastic variability 
has been studied by Lepine & Moffat (1999 ) in WR stars and by 
Markova et al. ([2005) in the Ha line of O stars. 

Clear but indirect evidence for hydrodynamic shocks in O- 
star winds is the observed X-ray emission and the superion- 
ization that the X-rays produce in the cool plasma (Cassinelli 
& Olson 1979). Theoretical work has attributed such shocks to 
the intrinsic "de-shadowing" instability of line-radiation driven 
winds (Owocki et al. [19881 Feldmeier [1995b . 

A first approximation to account for clumping in stellar wind 
models has been introduced about a decade ago (Hamann & 
Koesterke [19981 Hither & Millier [19991 Puis et al. 12006b . It is 
based on the assumption that the wind clumps are small com- 
pared to the mean free path of the photons. This "microclump- 
ing" approximation introduces a new parameter D to describe 
the factor by which the density in the clumps is enhanced com- 
pared to a homogeneous model with the same mass-loss rate. 
The most important consequence of accounting for wind clump- 
ing within the optically-thin clump limit is the reduction of the 
empirical mass-loss rates by a factor V75. This holds in princi- 
ple for all diagnostics that depend quadratically on the density, 
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thus including all emission lines in WR spectra, thermal radio 
emission, and the Ha line from O stars. However, the different 
diagnostics may be affected to a different degree, due to the pos- 
sible radial variation in the clumping factor D. 

A handle for estimating the clumping factor in Wolf-Rayet 
atmospheres is provided by strong emission lines with extended 
line wings, caused by the frequency redistribution of line pho- 
tons scattered by free electrons. Clumping reduces the rela- 
tive strength of these wings, compared to homogeneous mod- 
els (HiHier [T9911 >, because typical WR emission lines scale with 
the square of density, while electron scattering scales linearly. 
Detailed modeling of electron scattering wings therefore allows 
the density enhancement factor D to be estimated. 

Hamann & Koesterke (1998 ) typically found D-A for a few 
selected Galactic WN stars, corresponding to a volume filling 
factor of the clumps of fy = 0.25. More recent studies (Crowther 
et al. 120021 1 prefer even stronger clumping, with filling factors 
smaller than fy =0.1. 

For O-stars, Puis et al. ( 120061 1 conclude that the clumping 
factor D is about four times larger in the line-forming region, 
compared to the radio-emitting region far away from the star. 
Resonance lines, which are not affected by microclumping be- 
cause of their linear density-dependence, are mostly saturated 
in O stars and therefore not suitable for precise mass-loss de- 
terminations. However, Fullerton et al. ( 120061 ) exploit the fact 
that the far-UV spectral range observed with FUSE contains the 
P v resonance line, which is typically unsaturated in O-star spec- 
tra, and analyze this line with the SEI method (1987). They de- 
rived mass-loss rates that are about ten times lower than obtained 
from p 2 diagnostics with unclumped models. Very strong clump- 
ing and consequently low mass-loss rates were also obtained 
by Bouret et al. ( 120051 ) for two O stars analyzed with non-LTE 
model atmospheres. Such low mass-loss rates would have dra- 
matic consequences on stellar evolution and feedback. 

However, all these downward-revisions of mass-loss rates 
are based on the assumption that the clumps are always smaller 
than the mean free path of photons. For strong lines in dense 
winds, this approximation cannot be justified. In the present pa- 
per we relax this small-clump approximation and account for 
clumps of arbitrary optical depth. Significant work has already 
been done to unveil the effects of optically thick inhomogeneities 
in stellar atmospheres. 

Grey opacity in a stellar atmosphere where the photon mean 
free path does not exceed the scale of wind inhomogeneities 
was considered by Shaviv ( 1998). The term "porous atmo- 
sphere" was coined by Shaviv ( 120001 ) to describe a multi-phase 
medium that allows more radiation to escape while exerting a 
weaker average force. Such an atmosphere yields a consider- 
ably lower mass-loss rate for the same total luminosity and can 
explain the apparent existence of super-Eddington stars (Shaviv 
120001 ). Owocki et al. 020041 ) introduced a "porosity-length" for- 
malism to derive a simple scaling for the reduced effective 
opacity and used this to obtain an associated scaling for the 
continuum-driven mass-loss rate from stars that formally exceed 
the Eddington limit. 

Continuum opacity in inhomogeneous winds with optically 
thick clumps was studied by Feldmeier et al. ( 120031 ). They con- 
sidered X-ray lines emitted by an optically thin hot plasma that 
is attenuated in a fragmented cool stellar wind by strong contin- 
uum opacity from bound-free and K-shell photoionization pro- 
cesses. Oskinova et al. ( 120041 ) show that wind porosity can ex- 
plain X-rays line profiles and the weak wavelength-dependence 
of the effective opacity. Considering the specific case of spher- 
ical clumps, Owocki & Cohen (120061 1 argue that a large clump 



separation is required to fit the observed X-ray emission line pro- 
files. However, Oskinova et al. ( 120061 ) reproduced the observed 
emission line profiles in X-ray spectra of O-stars with a wind 
model assuming plausible model parameters. Addressing mass- 
loss rates inferred from the analysis of radio measurements, 
Brown et al. (120041) point out that optically-thick clumping leads 
to the reduction of multiple scattering and, consequently, photon 
momentum delivery. 

Line opacity in inhomogeneous stellar winds is modeled 
for the first time in the present paper. Considering the effect 
of clumping on UV resonance lines, Massa et al. ( 120031 ) and 
Fullerton et al. ( 12006b discuss how porosity can affect the for- 
mation of P Cygni lines and, in extreme cases, produce an ap- 
parently unsaturated profile for a line that would be extremely 
saturated if the wind material were distributed smoothly. 

In the following section (Sect. 2) we briefly review the line 
radiative transfer in the conventional optically-thin clump ap- 
proximation. In Sect. 3 we describe a statistical treatment for 
clumps of arbitrary optical thickness. In Sect. 4 both models are 
compared with observed mass-loss diagnostic lines in the spec- 
trum of ( Puppis. A brief parameter study is conducted in Sect. 5. 
Two representative examples of Wolf-Rayet model spectra are 
presented in Sect. 6, and a summary is given in Sect. 7. 

2. Radiative transfer in the limit of optically thin 
clumps 

Radiative transfer in an inhomogeneous medium is implemented 
in stellar-wind codes like CMFGEN and PoWR to a first approx- 
imation. This approximation holds in the limit that the size of 
the structures is small compared to the mean free path of pho- 
tons (microclumping). Moreover, it is assumed that the inter- 
clump medium is void. The density in the clumps is enhanced 
by a factor of D compared to the density p of a smooth model 
with the same mass-loss rate M\ hence, the volume filling factor 
is /v = D l , where D may depend on the location in the stellar 
wind, i.e. on the radial coordinate r. 

As there is no matter outside the clumps, rate equations are 
solved only for the clump medium, i.e. for the enhanced density 
pc = Dp. From the obtained population numbers, the non-LTE 
opacity and emissivity of the clump matter, Kc(Dp) and rjc(Dp), 
can be calculated. 

In the radiative transfer equation, the smooth-wind opacity 
and emissivity n{p) and rj{p) must be replaced for a clumped 
wind by 

Kf = fv K C (Dp) and n f = fy tj c (Dp) . (1) 

The factor fy accounts for the fact that only this fraction of a ray 
actually crosses clumps. This treatment of radiative transfer is 
sometimes termed the "filling factor approach" (Hillier & Millier 
[19991 . 

The atomic processes contributing to the opacity and emis- 
sivity scale with different powers of the density. For processes 
linear in density, fy and D cancel, but contributions scaling with 
the square of the density (bound-free emission, or free-free ab- 
sorption and emission) are effectively enhanced by a factor of D. 
(Of course, k and 77 are also indirectly affected by clumping via 
the population numbers, due to modified rates in the statistical 
equations.) 

Empirical mass-loss diagnostics are widely based on pro- 
cesses that scale with the square of the density. The wind emis- 
sion lines in Wolf-Rayet and O stars, including Ha, form in de- 
excitation cascades that are fed by radiative recombination. The 



L. M. Oskinova et al.: Neglecting the porosity of hot-star winds can lead to underestimating mass-loss rates 



3 



thermal radio emission from stellar winds is due to the free-free 
process. 

When the wind is clumped in the regions from which emis- 
sion emerges, the radiation flux is enhanced by a factor of D 
compared to a homogeneous model with the same mass-loss 
rate. Consequently, when a given (radio or line) emission is ana- 
lyzed with a model that accounts for microclumping, the derived 
mass-loss rate will be lower by a factor of Vd than from with a 
smooth-wind model. 

Puis et al. (12006b studied the O-type supergiant ( Puppis 
and found that the mass-loss rates derived from Ha and from 
IR/radio continuum emission can only be reconciled if the wind 
is strongly clumped (D * 5.5) already at a radius of 1.12/?* 
where the wind velocity is only aTOOkms -1 . 

Line absorption coefficients depend linearly on density. In 
the case of resonance lines, the same holds for the re-emitted 
("scattered") photons. Hence, P Cygni profiles are independent 
of microclumping. Of course, clumping might indirectly affect 
resonance lines via the ionization balance. In most cases of prac- 
tical importance, the resonance line belongs to the leading ion- 
ization stage, which is robust against such effects. 



3. Radiative transfer with clumps of arbitrary 
optical thickness 

The basic assumption of the "microclumping" approximation 
described in the previous section is that the clumps are small 
compared to the mean free path of photons, i.e. that they are op- 
tically thin. Given the large atomic opacity in the center of typ- 
ical spectral lines, this approximation is not generally justified. 
In the present paper we therefore relax this approximation in fa- 
vor of a more general treatment that accounts for the possibility 
of clumps being optically thick at some frequencies. We refer to 
this generalization as "macroclumping". 

3.1. Basic assumptions 

To introduce our treatment of porous stellar winds, we first con- 
sider clumps that are smaller than the photon mean free path at 
any frequency. Throughout the paper, we refer to these as "mi- 
croclumps". The non-LTE population numbers that determine 
opacity and emissivity are obtained by solving the rate equations 
for the clump matter density. The radiation field that enters the 
radiative transition rates is calculated with these opacities and 
emissivities, corrected for the clump filling factor (cf. Eq. (1)). 

As a next step, consider a few microclumps being merged to- 
gether to form one bigger "macroclump" without changing the 
matter density. The size of this macroclump may now exceed 
the photon mean free path at frequencies with large atomic cross 
sections. In other words, the macroclump can be optically thick 
in the cores of a few strong spectral lines, while it remains opti- 
cally thin at all other frequencies. 

The non-LTE population numbers depend on the population 
and de-population processes, i.e. on the collisional and radiative 
transition rates. Because we assume that the clump density does 
not change by assembling the macroclump, collisional rates are 
not affected and are the same as calculated in the microclump ap- 
proximation. The same holds for the radiative rates of optically 
thin transitions. 

The only transition rates that change in the macroclump are 
those radiative rates that are excited at frequencies where the 
macroclump is optically thick. The radiation field in optically 
thick lines may change drastically in the inner parts of a clump. 
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Fig. 1. Sketch of a clumped stellar wind. In a smooth wind, 
rays of a given observer's frame frequency encounter line opac- 
ity only close to the "constant radial velocity surface" (CRVS, 
thick shaded line). In a clumpy wind, assuming that the clumps 
move with the same velocity law as in the homogeneous wind, 
only those clumps interact with the ray that lie close to the cor- 
responding CRVS (dark-shaded circles). All other clumps are 
transparent (open circles) if the continuum opacity is small, so 
the wind is porous with respect to line absorption, even when the 
total volume is densely packed with clumps. 



However, in the optically thick limit, radiation is trapped inside 
the clump. Photons in these lines are immediately re-absorbed 
within the clump and do not escape to the interclump space. 

The energy of photons in optically thick lines must be redis- 
tributed elsewhere: the population of the upper line level will in- 
crease until it is de-populated by other (radiative or collisional) 
transitions. Hence, in principle the radiation field in any other, 
optically thin line may be modified. We neglect this second- 
order effect and focus on the macroclumping effect on strong 
lines. 

Photons in optically thick lines that emerge from the clump 
stem from its outer layers of optical depth about unity, irradi- 
ated from the interclump space. Therefore, even at these frequen- 
cies it is a reasonable approximation to adopt the radiation field 
from the microclumping model to evaluate the radiative transi- 
tion rates, i.e. to neglect the feedback of macroclumping effects. 

Summarizing these arguments, we find that 
f i) most of the transition rates are not affected at all by macro- 
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clumping (because they are collisional or excited by frequencies 

for which the clumps are optically thin), and 

( ii) for the few optically thick transitions, the average radiation 

field neglecting macroclumping provides a rough approximation 

to evaluate the radiative transition rates in the outer layers of the 

clumps. 

Therefore it is a reasonable first approach to use the non-LTE 
population numbers from the microclumping model and to ac- 
count for the macroclumping effects only in the "formal inte- 
gral". 

3.2. Macro-clumping approach 

For the radiative transfer in the "formal integral", we use a statis- 
tical approach that holds within the limit of large clump number. 
As in the microclumping case, we assume that the density inside 
the clumps is enhanced by a factor of D compared to a smooth 
model with the same mass-loss rate. The interclump medium is 
void, and the volume filling factor is /v = D . 

So far, no 3-D simulations of stellar winds are avail- 
able. Hydrodynamic 1-D simulations of wind instabilities by 
Feldmeier et al. (119971 ) indicate that the clumps take the form of 
shell fragments. This is supported by detailed modeling of ob- 
served X-ray emission line profiles (Oskinova et al. 120061 ). On 
the other hand, Dessart & Owocki ( 120051 1 suggest that clumps 
may have a similar lateral and radial scale. However, they em- 
phasize that the clump scale, compression level, and degree of 
anisotropy are still uncertain. 

To avoid the introduction of an excessive number of free 
parameters, we assume here that clumps are isotropic. They have 
a linear size / that is uniform at a given distance from the stel- 
lar center but may be a function of radius r. For simplicity we 
assume that a clump of volume P has diameter / for any cross- 
ing ray, thus neglecting any center-to-limb variation across the 
clump. 

It is important to realize that we are considering line ra- 
diation transfer in a moving medium, which is transparent in 
the continuum. The situation is illustrated in Fig.Q] A specific 
ray, when crossing the stellar wind, only encounters line opac- 
ity over a short part of its path, where the projected expan- 
sion velocity has just the right value to Doppler-shift the line- 
center frequency to the selected observer's frame frequency. In 
the framework of the Sobolev approximation, the observed line 
emission or absorption arises when the ray intersects the corre- 
sponding "constant radial velocity surface" (CRVS). In the case 
of a porous medium, the CRVS dissolves into patches. Adopting 
a clump/void density structure but an undisturbed velocity field, 
only those parts of the CRVS are effectively left that fall into the 
volume of a clump. Therefore, porosity effects can be especially 
important for line radiation transfer in expanding atmospheres. 

We assume that the clumps are statistically distributed, hav- 
ing an average separation L between their centers (L(r) varies 
with radial location); hence, the volume filling factor is fy = 
P/L?, or 



D 



From Eq. © it follows that 

«c = L~ 3 = D l r\ 



(2) 



(3) 



where «c denotes the number density of stochastically dis- 
tributed clumps. 

Lucy ( 120071 ) has shown that photospheric turbulence plays 
a decisive role in regulating the mass flux of a stationary flow. 



It is known from many studies that a high "microturbulence" is 
needed to reproduce the detailed shape of line profiles observed 
in O and WR star spectra. We identify this "microturbulence" 
vd with the typical velocity dispersion within a clump, caused 
by stochastic motions and the expansion of the clump itself. In 
our radiative transfer calculations, we therefore take the line ab- 
sorption coefficient with a Gaussian profile corresponding to the 
Doppler-broadening velocity l>d and evaluate the optical depth 
across one clump formally as for a static medium. 

We have already introduced the opacity kq within a clump, 
so we can write the optical depth across one clump with the help 
of Eq. (fl]i as 



T C = K C I = Kf Dl , 



(4) 



Alternatively, tc can be expressed by the density contrast D and 
the average clump separation L, 



(5) 



Note that for L — » the clumps become optically thin, approach- 
ing the microclumping approximation. 

The effective opacity of the clumpy medium is obtained 
in analogy to the usual opacity from atomic absorbers: 



K e ff = «c o" C 



(6) 



Here cr c is the effective cross section of a clump. "Effective" 
means that the geometrical cross section, I 2 , is multiplied by the 
fraction of photons that is absorbed when crossing the clump. As 
was derived in Feldmeier et al. (120031 1, 



(T C = I 2 (1 



) ■ 



(7) 



Using Eq. 0} for the number of clumps and Eq. (|7|) for the 
effective cross-section in Eq. (O, the effective opacity becomes 



*eff = (Diy' (1 - e- Tc ) . 
In the limit of opaque clumps, 
lim K eB = (Diy 1 . 



(8) 



(9) 



Noticing from Eq. (0]) that (Diy 1 = Kf/rc, we obtain a scal- 
ing of effective opacity with the opacity obtained in the micro- 
clumping approximation, 



1 -, 



KeB = K f 



Tc 



= Kf C n 



(10) 



The factor C macro thus describes how macroclumping changes 
the opacity, compared to the microclumping limit. Note that the 
small-clump approximation K e ff * Kf is recovered for optically 
thin clumps (tc <k 1). For optically thick clumps (tc~1), how- 
ever, the effective opacity is reduced by a factor C macro compared 
to the microclumping approximation. 

To understand the reduction of opacity in a medium with op- 
tically thick clumps, one can imagine that atomic absorbers are 
hidden in optically thick clumps and therefore not contributing 
to the effective opacity. Relatedly, Brown et al. (2004) noticed 
that this also leads to a reduction in radiative acceleration. 

While Feldmeier et al. (|2p03 ) used Eqs. ((6]) and (0 for their 
analysis of X-ray emission line profiles, Eq. (TTOb was derived by 
Owocki et al. (200 4) for the gray opacity. These authors define 
a "porosity length" h = L? /I 2 . Rewriting Eq. (TT0T > for the poros- 
ity length, one obtains K e ff = hr x (1 - <?~ Tc ). This can be com- 
pared to Eq. (A. 7) in Feldmeier et al. (120031 ) and Eqs. (29) and 
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(31) in Oskinova et al. ( 12004i i, where the opacity is allowed to 
be angular-dependent. 

We evaluate the correction factor C macro by specifying the 
clump separation L{r), which enters Eq. ((5]) in combination with 
the density contrast D{r). The reason for this parameterization is 
that D{r) must be specified for the density in the non-LTE rate 
equations. Note that the porosity length h - D 2 ^L alone is not 
adequate to define the model. 

Clumps are assumed to be preserved entities. This assump- 
tion is corroborated by the hydrodynamical models (Feldmeier 
et al. [19971 Runacres & Owocki l20U2l Dessart & Owocki l2005l 
that predict that the over-dense structures accelerate radially with 
u(r) and, in general, do not split or merge. As a result, the number 
density of clumps has to obey the equation of continuity, 



«c °c (r 2 v(rj) 1 



(11) 



Since the average clump separation is L — n c , we can write 



L(r) = L (r»W)" 3 , 



(12) 



where Lo is the free parameter of our macroclumping model. 
Note that r and Lo are in units of the stellar radius, and w(r) = 
v(r)/voo is the velocity in units of the terminal speed. 

Plausibly, Lo is close to unity. In the 1-D hydrodynamic 
model of Feldmeier et al. d 19971 1, the wind shocks evolve from 
photospheric seed perturbations that have periods close to the 
acoustic cutoff. As the latter is close to the wind-flow time R*/Voo 
for a typical O star, the predicted radial separation between the 
dense shells is approximately one stellar radius. 

When studying the effect of macroclumping on the line 
spectra, it is important not to spoil the reference continuum. 
Therefore we only allow for macroscopic clumping in the line- 
forming regions, but suppress it in the photosphere. Photospheric 
clumping would alter the radiation field in the line-forming re- 
gions (and thus the population numbers!) and would also change 
(enhance) the reference continuum for the line spectrum. 

Therefore we augment Eq. ( fT2l for the clump separation L 
by a factor s(tr oss ), which switches the macroclumping off at 
large Rosseland optical depth, with a linear transition between 
n and T2'. 



■S(TRoss) = 



TRo.ss)/(T2 - Ti) 



if TRoss 
if T"R oss 

else 



> 1"2 
< Ti 



(13) 



where we usually choose t\ — 0.3 and T2 = 1 .0. We have checked 
that the continuum flux (SED) is not modified by macroclump- 
ing effects with these parameters. 

We implemented the described modifications of the radia- 
tive transfer into the "formal integral" of the PoWR model at- 
mosphere code. At each integration step along each ray of given 
frequency, we calculated the clump optical depth tq from Eq. (O 
and evaluated the macroclumping correction factor (Eq. (fTOl)), 
which is applied to the opacity before the integration step is per- 
formed. Note that the same correction factor must also be ap- 
plied to the emissivity to preserve the non-LTE source function 
according to our assumption that the population numbers are not 
affected. 

3.3. Clump separation and the total number of clumps 

The equations in the previous section describe our approach to 
optically thick clumping. Only one additional free parameter, Lo, 
has been introduced compared to the microclumping approach 



(Eq.[T2l. For a specified stellar atmosphere model, the choice of 
Lq determines to what extend the clumps are optically thick at 
line frequencies. Within the limit of small Lo, the macroclump- 
ing effects become negligible. Realistic values of Lo can be es- 
timated from considering the total number of clumps that reside 
in the wind at any given moment. 

The number of clumps per unit volume, i.e. their number 
density, which is a statistical variable, follows from their average 
separation 

n c (r) = L- 3 = (l 3 r 2 w(rj) X . (14) 

Thus the total number of clumps that are found within the radial 
range from r\ to r2 is 

rn 4^- r>2 (j r 4^ 

Nc= n c (r) 4n r 2 dr = — _ = _(, 2 - , (15) 

J n ^0 J" W ^ L Q 

where ?2 - h is the flight time between r\ and r2 in units of 
the dynamical time scale, R„IVoc. For "beta velocity laws" v(r) = 
i>oo(l - l/^y 5 , analytical solutions for certain values of (3 are given 
in Hamann et al. (120011 1. For /3 = 1 it holds that 



t(r) 



-I 



dr 

w(r) 



ln(r - 1) + const 



(16) 



Let us consider, for instance, the range between r\ = 1.05, 
where the wind velocity is 5% of the terminal speed, and ri — 10 
stellar radii. The total number of clumps in this volume obtained 
from the previous two equations is 



N c = 178 L 



(17) 



Note that Nc scales with r-i more weakly than linearly; when 
restricting the volume to ri = 5, the number of enclosed clumps 
is still A^c = 105L 3 . 

Empirically, not much is known about the number of clumps 
in stellar winds. So far, only rough order-of-magnitude estimates 
can be made. Line profile variability may provide some relevant 
information. Lepine & Moffat (1 19991 1 monitored nine WR stars 
and found stochastic line profile variations in the form of narrow 
(«; lOOkm/s) emission sub-peaks on top of emission lines. They 
explained their data by 10 3 liVc~10 "blobs" being present at any 
time in the line emission region. Using Eq. ( TfTb this corresponds 
to Lo ~ 0.3. We will see in the next sections that a value of Lo 
between 0.2 and 0.5, implying a plausible number of 10 3 to 10 4 
clumps in the line-forming region, is high enough to produce 
significant macroclumping effects on the emergent spectra. 

Further support for selecting Lo in the range 0.1 to 1.0 comes 
from analizing X-ray emission lines. In Oskinova et al. (12006b . 
we reproduced the Chandra observations of O-star line pro- 
files by assuming that spherical shells are ejected once per dy- 
namic time scale fa = Rt/Va,, which corresponds to a radial 
separation of one stellar radius when v m is reached. Note, how- 
ever, that there are subtle differences between the "fragmented 
shell model" used for the X-ray modeling and the present paper, 
thereby preventing a direct comparison of the parameters. 

4. The impact of macroclumping on O star spectra: 
^Puppis 

As a representative test, we demonstrate the influence of macro- 
clumping on the spectrum of the O-type supergiant ^Puppis. 
First we calculate a non-LTE atmosphere model with the PoWR 
code (see Hamann et al. 2004 for details). Based on the model's 
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Fig.3. Effect of macroclumping on the Pv resonance doublet at 1 1 18/1 128 A for an O star. The spectrum of £Pup as observed 
by copernicus is shown for comparison (black, ragged line). The model parameters (see text) correspond to ^Puppis. We adopt 
a microturbulence velocity of 50km/s in the formal integral and do not apply rotational broadening; therefore the photospheric 
absorption components are not properly reproduced. Numerous weak spectral features in this range are due to iron. The usual 
microclumping modeling yields PCygni features that are too strong (blue, continuous line). With our macroclumping formalism, 
the line features are reduced to the observed strength (red, dotted curve). Macro-clumping was switched on at about the sonic point 
(lOkm/s) with a clump-separation parameter Lq = 0.2. 



population numbers, the emergent spectrum is then obtained by 
solving the "formal integral". This step is done twice, once with 
the usual code version that accounts for clumping only in the 
small-scale limit, and a second time with the macroclumping 
formalism as developed in the previous section. 

For the macroclumping simulation, we set the clump sepa- 
ration parameter to Lq = 0.2. According to Eq. ( fT71 ), this im- 
plies a 2.2 • 10 4 clumps in the line-forming region (within 10/?,). 
This number of clumps is large enough to justify our statistical 
treatment, and it compares well with the empirical estimates (cf. 
Sect.O- 

Macroclumping is assumed to start at about the sonic point 
(lOkm/s), together with microclumping. As model parameters 
for £Puppis, we adopt (compare Puis et al. 12006b logL/L = 
5.9," T» = 39 kK, logM /(A^yr 1 ) = -5.6. We use the com- 
monly accepted yS-law for the velocity field, but split it into two 
ranges with y8=0.9 from the photosphere to 0.6 v x and (3-A for 
the outer part of the wind. The terminal speed is 2250 km/s. 
The density enhancement factor is set to D — 10 in the super- 
sonic part of the atmosphere. For the chemical composition, 
we adopt 42/56/0.12/0.4/0.9/0.14 (mass fractions in %) for the 
H/He/C/N/O/Fe-group (Repolust et al. 2004). We also included 
phosphorus with solar abundance (6.15 X 10 -6 mass fraction). 

We concentrate here on two lines, Ha and the P v resonance 
doublet. The Ha emission is the most important diagnostic tool 
for O-star mass-loss rates. The Pv doublet, unlike other res- 
onance lines, is not saturated (because of the low phosphorus 
abundance) and therefore potentially useful for the empirical de- 
termination of mass-loss rates. FUSE observations of the P v res- 
onance doublet in the extreme ultraviolet caused the recent de- 
bate on a drastic downward revision of mass-loss rates (Fullerton 
et al. 120061 Bouret et al. l2005l 

The Ha emission line is strongly affected by microclump- 
ing, especially in the lower regions of the wind (see also Puis 
et al. 12006b . In contrast to Ha, resonance lines are not sensitive 
to microclumping because of their linear dependence on density. 



Moreover, P v is the leading ionization stage of phosphorus in 
the wind of typical O stars. This makes the line robust against 
details of the modeling, including the "superionization" by X- 
rays, which is an essential effect for the O vi and N v resonance 
lines. For the same reason, the P v resonance doublet is also in- 
sensitive to the second-order effects on the population numbers 
that we neglect in our macroclumping approach, as discussed in 
Sect.|3j] 

Figure[2]displays the synthetic spectra from both model ver- 
sions, without macroclumping use and with our macroclump- 
ing treatment, but otherwise identical parameters. The observa- 
tion of £ Pup is shown for comparison. The mass-loss rate of the 
model is tuned to the adopted density contrast D = 10 in order to 
obtain a good fit. We find that a model with the same mass-loss 
rate, but without microclumping, fails to reproduce the observed 
Ha profile by far (the line appears in absorption). 

Interestingly, the macroclumping effect is not significant for 
the Ha line. We tested different Lq and found no sensitivity for 
any realistic values. The obvious reason is that Ha is optically 
thin in the supersonic part of the wind from (, Pup, so accounting 
for macroclumping does not lead to a different empirical mass- 
loss rate from Ha fitting. 

When Bouret et al. d2005l) and Fullerton et al. d2006t re- 
cently analyzed the P v resonance doublet in O-star spectra, they 
encountered severe discrepancies in the mass-loss rates inferred 
from Ha. They suggested that the only way to reconcile these 
different diagnostics is to assume much higher clumping con- 
trast than assumed hitherto, implying a drastic reduction of the 
empirical mass-loss rates. 

This discordance is illustrated in Fig. [3] showing the EUV 
spectrum of f Pup (from the copernicus satellite) with the Pv 
doublet. The continuous line is the synthetic spectrum of the 
conventional model without macroclumping. The parameters are 
the same as in Fig. [2] i.e. perfectly able to reproduce the observed 
Ha emission. As can be seen, the predicted P-Cygni profiles 
of the P v doublet are much stronger than observed. To fit the 
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Fig. 6. Effect of macroclumping on the short-UV spectrum from a WN star of late subtype; the model parameters (see text) corre- 
spond to WR 78 (WN7h). The solid line is calculated with microclumping alone, and the dotted line is obtained with macroclumping 
in the formal integral (clump-separation parameter L = 0.5). 
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Fig. 2. Effect of macroclumping on the Ha- line from an O 
star. The observed spectrum of ^"Pup is shown for comparison 
(black, ragged line). The model parameters (see text) correspond 
to £Pup. The solid (blue) line results from a model that only 
accounts for small-scale clumping with a density contrast D = 
10, switched on at the sonic point (lOkm/s). The dotted (red) 
curve, partly identical with the solid (blue) curve, is obtained 
with macroclumping in the formal integral (also switched on at 
the sonic point with clump-separation parameter Lq = 0.2). The 
effect of macroclumping on the Ho- emission is small. 



1140 



Fig. 4. Synthetic Pv resonance doublet at 1 118/1 128 A for dif- 
ferent values of the clump separation parameter Lq (labels). The 
model parameters are for £ Pup (see Sect. 4), with vq = 50 km/s. 
A clump separation of zero recovers the microclumping limit. 
The larger the clump separation, the weaker the line becomes. 



lines in this spectral range are not affected by the macroclump- 
ing correction. When accounting for macroclumping, the Ha and 
P v lines in the observed spectrum of £ Pup can be fitted simul- 
taneously, using a "standard" clumping factor of D- 10 (filling 
factor fv — 0. 1). We conclude that accounting for macroclump- 
ing can result in a concordance of mass-loss estimates. 



Ha and P v lines consistently with the microclumping model, a 
higher density contrast D and a lower mass-loss rate would be 
required. 

Adopting the same model, we computed the stellar spec- 
trum with our macroclumping formalism (Fig.0. In the syn- 
thetic spectrum, the P v features are now drastically reduced, just 
to about the observed strength. Note that the numerous weak iron 



5. Influence of clump separation parameter, L , and 
velocity dispersion within a clump, e D 

As explained in Sect. 3.2, Lq is the essential free parameter of our 
macroclumping model. To demonstrate its influence, we con- 
sider the P v resonance doublet, keeping all other parameters as 
in the previous section. Blending iron lines are suppressed now 
for clarity. 



8 



L. M. Oskinova et al.: Neglecting the porosity of hot-star winds can lead to underestimating mass-loss rates 





!JL 



-I L. 



4000 



5000 



Fig. 7. Same as Fig. [6] but for the optical range 
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Fig. 5. Synthetic Pv resonance doublet at 11 18/1 128 A for dif- 
ferent values of the microturbulence velocity u D (labels). The 
clump separation parameter is Lo = 0.2 for all cases. All other 
model parameters are kept fixed. When the velocity dispersion 
across the clumps is decreased, the macroclumping effect be- 
comes more pronounced and leads to a weaker line profile. Note 
that the unrealistically strong photospheric absorption features 
become weaker with lower microturbulence velocities, as ex- 
pected from the curve-of-growth effect. 



We start by computing a model without macroclumping cor- 
rection, i.e. C m acio = 1 implying Lo = from Eq. (0. The result- 
ing profile is shown in Fig. [4] As the next step, we increase Lo 
stepwise while all other parameters are kept constant. Increasing 
Lo leads to an increased clump optical depth re- Therefore, the 
macroclumping correction factor C macro decreases, leading to a 
reduction of the effective opacity K e e, As can be seen in Fig.|4] 
higher values of Lo result in weaker lines. 

Significantly, the macroclumping effect depends not only on 
the number of clumps (defined by Lo), but also on the veloc- 
ity dispersion within the individual clumps. This is a specific 
consequence of the line radiative transfer in an expanding wind. 
Figure[5] shows model lines that are computed with a constant 
parameter Lq = 0.2, but for different velocities uq. As can be 



seen, the porosity effect of reduced line strengths is stronger for 
lower values vq. 

To understand this result, we considered the line absorption 
coefficient in a clump. With lower velocity dispersion within a 
clump, the Doppler broadening is smaller, and the line absorp- 
tion profile is narrower but peaks higher. Thus, the clump optical 
depth in the line core becomes larger, while it is smaller in the 
line wings. Higher optical depth in the line core leads to a re- 
duction of the effective opacity and a weakening of the line (see 
Eq. COll). 

This is in accordance with our interpretation of a porous con- 
stant velocity surface (see Sect. 3.2). As highlighted in Fig.Q] 
only clumps in the vicinity of the constant radial velocity sur- 
face contribute to the optical depth along the line-of-sight at a 
given frequency (see Sobolev 1 19471 1. In a wind with a mono- 
tonic velocity field u(r), the characteristic width of this vicinity 
is given by the "Sobolev length", V£,/(dv/dr) (e.g. Grinin 1200 11 1. 
For lower Doppler-broadening velocity uu, the clump is opti- 
cally thick for a narrower range of frequencies within the line. 
Consequently the number of clumps that can contribute to the 
line-of-sight optical depth is smaller, and the effect of porosity 
is greater. 

The velocity dispersion across a clump is well-constrained 
both from empirical and theoretical considerations. It has been 
known for decades (e.g. Ilamann 1980 Lamers et a!. 11987| i that 
detailed UV line fits for O stars require a "microturbulence ve- 
locity" of 50 to 100 km/s. The same holds for the emission line 
spectra of Wolf-Rayet stars. In our picture of a clumped wind, 
this small-scale Doppler broadening is identified with the typi- 
cal velocity dispersion within one clump. Moreover, it has been 
noted repeatedly from monitoring line-profile variability that 
narrow features that appear transiently are never narrower than 
50 to 100 km/s (e.g. Lepine & Moffat [19991 . 

Clumps can only survive over the dynamic flow time when 
their internal velocity dispersion is small compared to the wind 
velocity, <k v m . Close inspection of the 1-D time-dependent 
model by Feldmeier et al. i 19971 reveals that, within the den- 
sity enhancements, the velocity does not differ by more than 
w 100 km/s. Since uu is constrained to 50 to 100 km/s for O and 
WR stars, the clump separation Lo is the basic free parameter of 
the macroclumping model. 
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Fig. 8. Effect of macroclumping on the UV spectrum from a WN star of early subtype WR 7 (WN4-s) with parameters as given in 
the text. The solid line is calculated with microclumping alone, and the dotted line is obtained with macroclumping in the formal 
integral (clump-separation parameter Lq = 0.5). 
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Fig. 9. Same as Fig. [8] but for the optical range. 



6. Effect of macroclumping on Wolf-Rayet spectra 

In this section we test the effect of macroclumping on the emer- 
gent spectra for representative models of a WNL star (i.e. a late- 
subtype Wolf-Rayet star of the nitrogen sequence), and a WNE 
(i.e. early-subtype Wolf-Rayet) star with strong emission lines. 

The two examples are selected from the PoWR grid of mod- 
els for WN stars. First we take a model with 20% hydrogen 
("WNL grid" model WNL07- 1 1 , Hamann et al. [2004b with para- 
meters corresponding to a late subtype (WN7): T„ = 50 kK, 
logR t /R = 1.0 (see Hamann et al. 120041 for the definition and 
role of this "transformed radius"). We use the standard /J-law for 
the velocity field with [5-1 and Voo = lOOOkm/s. The microtur- 
bulence velocity for Wolf-Rayet models is set to u D = lOOkm/s. 
In contrast to the published grid (www.astro.physik.uni- 
potsdam.de/~wrh/PoWR/powrgridl.html), we choose here a 
density contrast of D = 10. The adopted composition is 
78/20/1.5/0.01/0.14 for the mass fractions (in %) of the 
He/H/N/C/Fe-group. Except for the different (but scalable) lumi- 
nosity, these parameters correspond to the Galactic Wolf-Rayet 



star WR78. Compared to the previous example (£Pup), this 
model is slightly hotter and has a denser wind. In this example 
we have chosen a clump separation parameter of Lq = 0.5, i.e. 
the total number of clumps within 10 R t according to Eq. ( fT71 ) is 
1424. 

The UV spectrum (Fig. [6]) shows numerous iron lines, partly 
in absorption and partly in emission at longer wavelengths. 
There are also a number of strong emission features, mostly as 
part of the P Cygni profiles. The highest emission peak is due to 
the He n 1640 A transition. Comparing again the emergent spec- 
trum with and without the macroclumping effect, we notice that 
generally all features, emission and absorption, are reduced in 
strength. Only the weakest lines are not affected, whereas the 
strongest emissions are cut down by more than 50%. 

In the optical range similar reductions are found for strong 
lines, but Fig.|7]also demonstrates that the effect differs between 
individual lines of similar flux. Obviously, to restore a spectral 
fit with the macroclumping model, one would have to increase 
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the mass-loss rate and possibly to re-adjust other parameters and 
abundances. 

Our example with the densest wind corresponds to an early- 
subtype WN star with strong lines, WR7 (WN4-s). The para- 
meters are logL/L = 5.3, T* = 112kK, logR t /R G = 0.3, no 
hydrogen, v x = 1600 km/s, D=10. The chemical composition is 
the same as for the WNL model from Sect.|6l but without hydro- 
gen. Except for the higher value of D, this model corresponds to 
the PoWR grid model WNE14-18 (Hamann et al. [2004T >. For the 
macroclumping correction, we again choose a clump separation 
ofL = 0.5. 

The emergent spectra including macroclumping effects are 
shown in (Figs. [8] and|9]l. As expected, the macroclumping re- 
duction is most pronounced for the strongest emission lines. As 
in the O star example, it is apparent that higher mass-loss rates 
must be chosen in order to maintain the lines as strong as they 
have been without macroclumping correction. 

Interestingly, even the "iron forest" is now significantly af- 
fected by the porosity effect. This is because many of the iron 
lines in dense WR winds are quite strong (peak intensities up to 
four times the continuum, cf. Fig.|7]i, and form throughout the 
wind like any other lines of comparable strength. 

Figure[lO]demonstrates the influence of the clump separation 
parameter Lq on the prominent Hen emission line at 4686 A. 
The emission decreases with increasing Lq. This porosity effect 
is already significant for relatively small clump separation, like 
L Q = 0.1. 

In the process of a spectral analysis, one would have to 
compensate for this macroclumping effect by adopting a higher 
mass-loss rate. Changing M, however, shifts the ionization bal- 
ances. Given the different sensitivities of individual lines on the 
macroclumping effect (cf. Figs. 6-9), element abundances may 
become involved as well. It should also be noticed that the den- 
sity contrast D(r) may actually vary with radius. Thus it is by 
no means straightforward to achieve a new consistent spectral fit 
for WR stars with macroclumping, so we must leave this task for 
future work. 



7. Summary 

We have incorporated the effect of wind porosity into the "formal 
integral" of the PoWR non-LTE atmosphere code. In a statistical 
treatment, optically thick clumps lead to an effective reduction of 
the opacity. The feedback of macroclumping on the population 
numbers is neglected as a second-order effect, especially in the 
case of resonance lines from a leading ionization stage. 

The effect of porosity on line formation differs drastically 
from that on the medium with continuous opacity. The effect of 
porosity on lines is enhanced due to the small extent of the zone 
where a ray is in resonance with the line frequency ("scatter- 
ing zone", "constant radial velocity surface"), while outside this 
zone the line opacity vanishes anyway due to the Doppler shift. 
Although a large number of clumps may be present in the wind, 
only a limited number of these clumps contribute to the line 
opacity along a specific ray. The statistical treatment of poros- 
ity is corroborated by the time-averaging nature of observations. 
Even when the number of clumps that are present in the scat- 
tering zone at a given moment of time is small, there is a large 
number of clumps passing through this zone during a typical in- 
tegration time of the observation. 

Following the predictions of time-dependent hydrodynamic 
models, we assume that clumps are preserved entities and that 
there is no velocity gradient inside the clumps. The velocity dis- 
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Fig. 10. Synthetic He n line at /14686A in dependence on the 
clump separation parameter Lq (see labels), uo = 50 km/s for 
all lines. A clump separation of zero recovers the microclump- 
ing limit. The larger is the clump separation, the weaker the line 
becomes. 



persion within individual clumps is described by a "microturbu- 
lence" velocity that broadens the line opacity profile. For sim- 
plicity, the clumps are assumed to be isotropic. 

We modeled the emergent spectrum of the O-type supergiant 
^Puppis and demonstrated the influence of macroclumping on 
the mass-loss diagnostic lines Ha and P v. Whereas the optically 
thin Ha line is not affected by wind porosity, the P v resonance 
doublet becomes significantly weaker when macroclumping is 
accounted for. Hence, neglecting macroclumping can lead to un- 
derestimating empirical mass-loss rates. In the case of ^Pup, 
the discrepancy between the Ha and Pv diagnostic found by 
Fullerton et al. (120061) can be entirely resolved with the macro- 
clumping modeling, without a downward revision of the mass- 
loss rate. 

We also studied the influence of macroclumping on the emer- 
gent spectra for two representative Wolf-Rayet star models, cor- 
responding to a WN star of late subtype and a WN star of early 
subtype with strong lines. The results show that basically all 
spectral features, in both emission and absorption, become sig- 
nificantly weaker when macroclumping is taken into account. 
Strong lines are typically reduced by a factor of two in intensity, 
while weak lines remain unchanged. 

In spectral analysis, this weakening of lines must be com- 
pensated for by adopting a higher mass-loss rate. For emission 
lines that scale with p 2 because they are fed by recombination, 
macroclumping is counteracting microclumping in its effect on 
the empirical mass-loss rates. 

Most important, wind porosity affects the P-Cygni profiles 
from resonance lines, in contrast to microclumping, which does 
not influence opacities that depend linearly on density. This can 
resolve the reported discrepancies between resonance-line and 
recombination-line diagnostics. Future work is needed to quan- 
tify this effect for individual stars. 
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